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SUMMARY 


An interactive digital computer program for modal analysis and gain estimation 
for eigensystem synthesis has been written. Both the mathematical and the opera- 
tional considerations are described; however, the mathematical presentation is 
limited to those concepts essential to the operational capability of the program. 

The program is capable of both modal and spectral synthesis of multi -input control 
systems. It is also user friendly with scratchpad capability and dynamic memory and 
can be used to design either state or output feedback systems. 


INTRODUCTION 

Srinathkumar (refs. 1, 2, and 3) developed several algorithms for the synthesis 
of control systems by using the concept of eigensystem (or pole) placement. He 
later programmed one of the algorithms (modal synthesis) for use on the Control Data 
Corporation (CDC) CYBER computer system and solved several example problems. But 
that program was written as a subroutine and was intended to operate in a "batch" 
mode. Moreover, the program was not formally documented for public use and did not 
include the algorithm for spectral synthesis. 

This document describes a digital computer program which is based on the work 
of Srinathkumar. This program takes advantage of the interactive capability of an 
on-line terminal within the modern computer complex. It also includes those algo- 
rithms for both modal and spectral synthesis as well as numerous other options. The 
program is called IMAGES, which is an acronym for Interactive _Modal Analysis and 
Gain ^Estimation System. 

This paper describes the interactive program and the mathematical basis for the 
major algorithms. The reader is referred to references 1 and 2 for those details 
omitted here. The present paper is intended to introduce the interactive program 
IMAGES to the scientific community and to serve as a user's guide. Hence, program 
structure, interactive discourse, and other user-related information are emphasized. 
An example problem is included to clarify operational procedures and to demonstrate 
ease of application. It is also intended to clarify certain points not clearly evi- 
dent within references 1 and 2. 


SYMBOLS 

A system matrix (numerical subscripts indicate particular partitions of this 

matrix) 

B input matrix (numerical subscripts indicate particular partitions of this 

matrix) 

C matrix defined in equations (23) 

C i'Ck matrices defined by equations (10b) and (15b), respectively 
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matrix defined in equations (23) 

vector with the rth entry equal to unity and the remaining elements equal 
to zero 

matrix defined in equations (9) 
vector defined in equation ( 21 ) 
matrix defined in equations (9) 
vector defined in equation (21) 

vector defined in equation (21) 
order identity matrix 
indices 

feedback matrices 

matrix defined in equation (16) 

number of inputs 

vectors defined in equation (16) 
matrix of (k - 1 ) eigenvectors 
number of states 

matrix defined below equations (25) 
matrix defined by equation (18) 
matrix defined by equation (14) 

number of outputs in output feedback vector 

• . . (k ) 

implies r is an element of the set A 

matrix defined in equations (9) 

time 

input vector defined by equation (2) 
u is an m x 1 real vector 
external reference input 
modal matrix (matrix of eigenvectors) 
eigenvectors 
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w. ,w vectors representing partitions of v. and v, , respectively 

i K X K 

x state vector defined by equation (1) 

n 

x G R x is an n x i real vector 

y output vector defined by equations (23) 

z.,z vectors representing partitions of v. and v, , respectively 

1 X XX 

z, c z is an m x 1 real vector 

l 

(k ) 

A set of indices defined in equation (17) 

Sz perturbed quantities 

K 

A eigenvalue matrix 

A , A eigenvalues 

l k 

a, a scalar quantity indicating limiting value 

K 

Superscripts : 

(1),(2),...,(k) indicates an iterative step 
T matrix transpose 

-1 matrix inverse 

Subscripts: 

o initial value 

ref reference value 

Uppercase letters denote matrices; matrix subscripts indicate partitioned quan- 
tities. A dot over a symbol denotes a derivative with respect to time; det is an 
abbreviation for determinant. A tilde (~), caret (*), or bar ( - ) over a variable or 
matrix denotes a transformed quantity. 

GENERAL DISCUSSION 
Modal Synthesis 

Srinathkumar (refs. 1, 2, and 3) both developed and programmed an algorithm for 
the synthesis of control systems by using the concept of eigensystem (eigenvalue/ 
eigenvector) or pole placement. He derived constraining relations for calculating 
those feedback gains which will provide the desired modal characteristics of a linear 
system with multi-inputs. Although this analysis is described in detail in refer- 
ence 1 and to a lesser degree in reference 2, multivariable synthesis by eigensystem 
assignment is briefly discussed herein. 
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The analysis begins with the controllable system 


x = Ax + Bu 


(1 ) 


where x £ R n and u £ R 1 ". "then, if B is assumed to be full rank, the state 
variable feedback law can take the form 


u = Kx + u 


( 2 ) 


where K is the feedback matrix, and u p is an external reference input. Hence, 
the problem is to select a K so that the closed-loop system matrix A + BK 
satisfies 


[A + BK ] v . = X.v. 

1 li 


(i — 1 , 2, » « . , n) 


(3) 


where X^ is the ith eigenvalue, and v^ is the corresponding eigenvector. 

Now, since m represents the number of inputs, equation (3) can be partitioned 
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where A-j-j, B^, and K-j are m x m matrices, and the remaining matrices are com- 
patibly dimensioned. Moreover, £ R m , and 


t r T i T 

/. = z. i w. 

1 Li i ij 


where is the design mode shape (user-specified part of mode shape). Then, com- 

pleting the multiplication of these partitioned quantities yields 


[A 11 

+ B K lz. + 
1 1 1 

[A 12 + 

B K ] w . = X . z . 
1 2 1 11 

(5) 

[A 21 

+ B K ] Z . + 
2 1i 

[A 22 + 

B K ]w. = X.w, 
2 2 1 11 

(6) 


Hence, it follows from equation (5) that 


K z. + K w. = B 1 [ X . z . - A z. - A w. ] 
1 l 2 i 1 i i 11 i 12 i 


(7a) 
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(7b) 


Vi + Vi * *2 


1 Vi - A 22“i 


- W 


provided B.j is either nonsingular or can be made nonsingular by reordering the 
state variables. That is, on combining equations (7a) and (7b) to eliminate K-) 
and K. 2 , 

[\.I - F]w. = [G + \.S]z. (i = 1, 2, n) (8) 
x n-m x xx 


where I n _ m is an (n - m)th-order identity matrix, and 


S 


B 2 B 1 1 




o 

II 

> 

- SA„ 

21 

1 1 

II 

> 

- SA 

22 

1 2 


(9) 


Equation (8) is the fundamental eigenvector-constraining relationship used by 
Srinathkumar in his algorithm for the synthesis of multi -input control systems. It 
constitutes a set of n-m linear equations in n unknowns, each of which repre- 
sents an element of an eigenvector. Then, if is not an eigenvalue of P, the 

n-m elements of w^ are uniquely related to the remaining m eigenvector ele- 
ments, each of which may be assigned arbitrarily. These assigned elements correspond 
to the z^ vector, while the n-m elements which correspond to the w^ vector 
must be computed. That is, 


w = C z (i = 1 , 2, ..., n) (10a) 

x i x 


where 


C. 

x 


[VI 
x n- 


■m 



[G + VS] 


(10b) 


is defined as that modal coupling matrix corresponding to Hence, after all 

n eigenvalues are assigned, up to nm eigenvector elements can also be arbitrarily 
assigned through the use of state variable feedback, provided the resulting modal 
matrix 



is nonsingular. The reader, however, must be aware that those choices of and 

zi depend on the desired characteristics of the control system. 

Equation (8) can be applied in a variety of ways. Aside from that option of 
synthesizing those eigenvectors achievable by state feedback, the analysis can be 
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extended to the problem of output feedback. Moreover, equation (8) can be used to 
evaluate the achievability of certain preconceived objectives. For example, there 
are occasions where those eigenvalues and eigenvectors chosen to achieve certain 
design objectives will not satisfy those constraints imposed by equation (8) on the 
choices of and z^. In these instances, it will be necessary to perturb 

slightly either Aj_ or z^ to achieve these design objectives. The option which 
provides this capability is the spectral synthesis option. Equation (8) can also be 
used to establish the parametric equivalence between the nonunique feedback matrix 
K and the arbitrary modal entries. Each of these options is included within IMAGES 
and is described briefly in this paper. But, the reader is, once again, referred to 
references 1 and 2 for additional information. 

The feedback matrix K is easily evaluated by returning to equation (3) 
or (4). If a nonsingular modal matrix V and an eigenvalue matrix A should be 
chosen to satisfy equation (8), then A.: will not be an eigenvalue of F, and the 

A -*• 

closed-loop system matrix A = A + BK is uniquely determined by 


A = 


A . A 


11 ! A 12 


A 

*A 

L A 2l! A 22 


= VAV 


-1 


( 11 ) 


Moreover, that feedback matrix K required to yield this closed-loop matrix A is 
computed using the relations 


K 1 “ B l 1 [ A 1 1 ~ A l 1 ] 


(1 2 ) 


k 2 “ B 1 [ A 1 2 — A 1 2 1 


where 


k = [ K) : k 2 ] 


Hence, given a modal matrix V and an eigenvalue matrix A, it is possible to 
evaluate the feedback matrix K, provided the matrices A and B are known. 


Spectral Synthesis 

The spectral synthesis option is more complicated. This option also assumes 
that one may be given those matrices A and B along with the desired eigenvalues 
and the eigenvector elements z^. In this instance, however, it is required to 
synthesize a nonsingular modal matrix of eigenvectors, each of which not only satis- 
fies equation (8) but also is mutually independent of the others. The latter condi- 
tion is satisfied provided the kth eigenvector v^ is chosen so that 


(k-D 
R v. 


(13) 
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where rOc- 1) i s the projector of that subspace spanning the previous (k - 1 ) eigen- 
vectors. That is, if 


N 





represents the (k - 1 ) eigenvectors, then 


R a-n 


T 

= N[N N] 


-1 


T 

N 


(14) 


represents the projector of that subspace spanning these eigenvectors, and equa- 
tion (13) ensures a nonsingular modal matrix. 

It is also essential that each eigenvalue \ k not be an eigenvalue of F if 
equation (8) is to be satisfied. If the desired eigenvalue should not combine with 
F so that 


det [ \ I - F] # 0 

k n-m 


(15a) 


then 


°k * [ V„-m - IG + V 1 ( ' 5b> 

is not defined. Spectral synthesis is an iterative procedure which avoids this 
difficulty by use of equation (8) to guarantee a solution for w^ which will combine 
with z k so that v k satisfies equation (13) and k k satisfies equation (15a). It 
requires that each eigenvector be evaluated sequentially and allows for perturbations 
within both ^ and z k . 

The algorithm depends on constructing an n x n matrix M^) such that 
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In equation (16), both m^ and mP^ are vectors which must be evaluated, and 

cl D / v \ , _ , 

the superscript (k) refers to the kth index chosen from a subset A v ' of indices 
1, 2, ..., n. These indices represent the corresponding values of and z^, and 

the superscript represents the kth permutation of these indices. Hence, a particular 
index will no longer be resident within the subset A^) once it is used in the con- 
struction of M^ 1 ), m( 2 ), ..., M . 
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Now , M< k ) is defined so that equation (13) may be satisfied without evaluating 
R explicitly . The construction is such that v^ k > is transformed into canoni- 
cal form by 


(k) < k > 
M v. 
k 


Vr 


(r C A 



(17) 


where * 0, and e r is an n vector with the rth entry equal to unity and all 

other entries equal to zero. Then, both the vectors m^ k * and m£ k * must be con- 
structed so that C1 L 


°k = 



Hence, r is both the rth element of v^ k) and the rth row of M (k ^. With these 
definitions, an n x m Q matrix is defined as 


0 (k_1) - M (k “ 1) M (k - 2) M {1) 

Q = M M . . . M 


(k > 1) 


(18) 


where 


Q 


( 0 ) 


= I 

n 


(19) 


Hence, 



n (k-D 

Q 


( 20 ) 


The construction of M^ k ^ and hence Q^ k ^ is such that a nonsingular modal 
matrix is assured. The definition of Q( k_1 ) also permits considerable freedom in 
the choice of z k and \ k from the list of desired values denoted by and X- . 

These values may be permuted at will to find the best combination of independent 
eigenvectors to achieve the desired design objectives. The procedure also permits 
perturbations within both z k and \ k to ensure that both equation (13) and equa- 
tion (15a) are satisfied. These and other points should become clearer as the 
remaining steps are outlined. 

The v k th eigenvector is synthesized by using the rth row of equation (18). It 
is first required to compute 


[g 


(k) iT 
r -* 


f (k-1) 
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+ h (k_1) C 
r k 


where 


f (k-1) 
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1 (k-1 ) 

' h 
i r 
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( 21 ) 
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is the rth row of the transformation Q^ k 1 ^ , and h^ -1 ) is an (n - m) row vector, 
while is evaluated using equation (15b). Then, 



(k) -iT 

J z 

r J 


k 


( 22 ) 


which must not be equal to zero. (See eq. (17).) In the event 0^=0, the user 
must either select another r, perturb z^ to z^ + 6z^, or return to perturb X^ 
once again. (See eqs . (15).) Normally, the order of application is as listed (i.e., 
select another r, etc.), with all remaining values of r being tried before per- 
turbing either z^ or X^. Hence, on finding those values of and X^ which 

will satisfy those conditions specified, w k can be evaluated by using equation (8). 
The results will be a v^ which is independent of those eigenvectors computed 
earlier. Moreover, subsequent eigenvectors computed in the same manner will also be 
independent. But, must be constructed from equation (17) before these steps 

are repeated for subsequent eigenvectors and eigenvalues. 


Output Feedback Control System 

The synthesis can be extended to an output feedback system. For example, con- 
sider systems of the form 


x = Ax + Bu 
y = Cx + Du 


(23) 


where y is an r x 1 output vector with r < n. Then, an approximation to the state 
feedback law (eq. (2) ) can be found using the control law 


u = Ky + u p 


(24) 


That is, on combining equations (23) and (24), the closed-loop system takes the form 


x = [A + BK]x + Bu 


y = [C + DK]x + Du 


(25) 


where K = PKC, B = BP, D = DP, and P = [l m - KD] -1 . Hence, P is the equivalent 
feedforward matrix, and K is the equivalent state feedback matrix. Then, for every 
state feedback matrix K derived by using the control law u = Kx + u , there exists 
an output feedback approximation 


~ _ + 
K = K[C + DK] 


( 26 ) 


where the superscript dagger (t) indicates the generalized inverse. 
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These algorithms are more general than is apparent from the preceding descrip- 
tions. Throughout the preceding discussion, the emphasis has been on real variables. 
It is, however, possible to extend the theory to complex variables. In fact, 
Srinathkumar has shown in references 1 and 2 that complex pairs can be treated using 
the same equations with very little modification. Those modifications have been 
included within the program and will not be discussed herein. Both real and complex 
eigenvalues and eigenvectors may be entered into the program without regard for 
order. The user, however, must be aware that the expected number of eigenvalues is 
diminished by one for each complex number entered. This information will be espe- 
cially useful when running IMAGES. 


PROGRAM DESCRIPTION 

The digital computer program IMAGES takes advantage of the interactive capabil- 
ity of the modern computer complex in the modal analysis and gain estimation for 
control system synthesis. The emphasis is on simplicity with a minimum of demand on 
the user. "Control" variables and other decision-making logic have been reduced to 
a series of numbered lists from which the user must select one. Moreover, dynamic 
memory is utilized to eliminate any need to redimension any DIMENSION statement. 

These and other ideas have been combined with those algorithms of reference 1 as 
described herein to develop a truly interactive program with scratchpad capability. 

The program was developed using that MODAL synthesis package written by 
Srinathkumar. That program, however, suffered greatly in two important points. 

First, it was based only on equations (1) through (12) and equations (23) 
through (25). This limitation eliminated spectral synthesis as a viable option 
within the program. As an interactive tool, this option may not only be useful 
but also essential. Second, MODAL was a card-oriented program (subroutine) 
possessing numerous options which required the user to resubmit the job for pro- 
cessing with each change in either options or model. The interactive program 
eliminates these time-consuming requirements. 

Srinathkumar' s program (or, more accurately, subroutine) was extremely flexible 
and well written. It was capable of estimating either eigenvectors, eigenvalues, 
or gains for either state feedback or output feedback control systems given those 
matrices indicated by either control system. (See eqs. (1) and (2) and eqs. (23) 
and (24).) Although the modal matrix would normally be synthesized by using equa- 
tion (10), another option computes this matrix, along with its associated matrix of 
eigenvalues A, for the given system matrix A. The gain matrix K was estimated 
by using equation (11). It was also possible to compute K on entering both the 
modal matrix and its associated A matrix as input to the program; the gains would 
be estimated using both equations (10) and (11). Other options were also available 
to the user. But, the semimodular nature of the batch -oriented subroutine MODAL 
imposed a requirement by the control logic for several different codes, each of which 
could assume several different values. 

Those code words used in MODAL along with their possible values are given in 
table I. As can be seen, one can be easily confused by the many possible combina- 
tions required to solve a particular problem. These codes, however, were necessary 
if the program was to be both efficient and useful in a batch computational mode. 

They are also essential to a lesser degree in the interactive program. They acti- 
vate that control logic essential to provide the flexibility needed to accomplish 
many different objectives by using the same source code. 



IMAGES uses a modified version of Srinathkumar' s program. The interactive pro- 
gram does not require the user to provide that FORTRAN logic needed for a calling 
program. Neither is it required to maintain a close watch over returning output as 
well as an ever-changing data deck. Punched card input has been replaced by either 
keyboard or disk file entries. Control words and other variables associated with 
control logic are automatically assigned within IMAGES and are no longer the respon- 
sibility of the user. These objectives are accomplished through a series of question 
and answer diagnostics, which are described in greater detail later. These questions 
and answers determine which codes will be set along with the appropriate entry point 
within MODAL. This has been achieved by restructuring the existing subroutine by 
using entry points. Hence, the original subroutine MODAL has essentially been decom- 
posed into four modules: (1) an input module, (2) a synthesis module, (3) a gain 

module, and (4) a response module. These modules are described graphically by the 
flow charts shown in figures 1 through 4. They represent the major entry points 
within the present version of MODAL. 

These charts are somewhat self-explanatory. The various entries and operations 
have been grouped into those modules which best describe their functions. Addresses 
of selected READ'S and other transfer points have also been included to aid the 
user who may have a listing of the program. (Both the program and the listing are 
available through COSMIC - Computer Software Management and Information Center, 

Suite 112, Barrow Hall, University of Georgia, Athens, GA 30602.) For example, the 
numbers 1030, 1040, 1037, and 1045 in the gain module (see fig. 3) indicate addresses 
of transfer points within MODAL. Numbers less than 1 0 on the figures are associated 
with those codes selected automatically by the various interactive responses. These 
responses are shown in figure 5 and are discussed later. Hence, that path taken 
within each module depends on the particular option chosen by the user. 

One should not be confused by the subroutine names but should be continuously 
aware of the differences between modal and spectral synthesis. The spectral synthe- 
sis subroutine SPECTRL - which is called from the synthesis module within MODAL - was 
not included within the original version of the pole placement package. It is based 
on that discussion centered around equations (13) through (22) and depends on subrou- 
tine SYNTHS. The original subroutine was named to reflect its primary function of 
modal synthesis. But, with this addition of the spectral synthesis algorithm, both 
the spectral synthesis and the modal synthesis are options within the synthesis 
module. Hence, MODAL is now an executive subroutine, and SYNTHS is that subroutine 
on which both algorithms depend. Both algorithms require the user to provide m 
elements of each eigenvector along with its corresponding eigenvalue. Modal synthe- 
sis will not adjust these design values. But, spectral synthesis permits the user to 
perturb these values in such a way as to satisfy within a certain tolerance (denoted 
by a) the condition of mutual independence. 

The remaining subroutines will neither be described graphically nor discus- 
sed herein, since the program listing is not included as a part of this document. 
These subroutines are, however, listed in table II along with a brief description 
of their functions. Hence, the remaining discussion is restricted to operational 
considerations. 

IMAGES is a user-oriented, interactive, digital computer program. Many differ- 
ent questions, messages, and other diagnostics are included within this program. The 
various matrices and other critical variables are printed on the terminal screen as 
they are generated, with sufficient pauses to permit the user to examine his or her 
results. In many instances, the user is permitted to modify these values. This is 



particularly true with those matrices and vectors serving as input. The user may 
also permute the A, B, C, and D matrices. These, as well as the several other 
options, provide considerable flexibility in the design and analysis of control 
systems . 

The interactive discourse is given in the flow chart shown in figure 5. Most of 
these diagnostics are included within the executive program IMAGES and the executive 
subroutine MODAL. Other diagnostics are also present, but these are not considered 
to be essential to the discussion. That path taken by the program in response to a 
given command is indicated by its number. Clearly, the user has many different 
choices. 

The program begins with the "master menu." It is indicated by the bold border 
in figure 5. All other blocks, with only one exception, are subordinate to this 
block and must ultimately return to it. Hence, each of those options within the 
master menu is a major program branch. Moreover, options within the submenus on one 
or more of those program modules are described elsewhere herein. 

The user is essentially free to choose among the several options as needed. 

Each is, hopefully, self-explanatory, and the order of selection is usually unimpor- 
tant. But, the user should be aware that the printed order of options is suggestive 
of that sequence of decisions a researcher would ordinarily make during the synthesis 
of an eigensystem for gain estimation. This is particularly true for new problems. 
The user must always "begin problem," "enter input," "enter computational mode," and 
"stop," or terminate, the program. The remaining options are utilities which enhance 
the overall flexibility and operational capability of the program. The "save prob- 
lem" option, in particular, allows the user to automatically save all those program 
variables needed to restart the program in the solution for a particular problem at 
some future date without penalty. This is provided as a safety feature as well as 
for computational flexibility. The save option automatically files all the current 
program variables and other data onto STRTUP - a random-access mass-storage disk 
file - under a name chosen by the user for future reference. It is also provided 
within each of the two remaining branches: (1) enter input and (2) enter computa- 

tional mode. 

On beginning a problem, the user must first specify whether the problem is an 
old or a new problem and supply a file name, which can be as many as seven charac- 
ters. (Note: A first-time user must first establish STRTUP as a viable storage file 

within the computer complex.) The program will then determine if the name has been 
used previously. If so, the user must decide if the name is to be retained. In the 
event his problem is to be a new example, the user may either purge the existing file 
or rename the new file. If the problem is an old one, the program will automatically 
find and load all those values associated with all steps completed up to the most 
recent save command using that file name. 

System matrices can be reviewed, modified, or both reviewed and modified on 
entering the enter input module. The user may also choose to go directly to the 
computational mode. If the problem is a new one, however, the user must enter those 
matrices required for his problem. Diagnostics specify the expected format for these 
new data. Provision is also made for correcting the data. All these diagnostics are 
self-explanatory and follow in a natural sequence to minimize the demand on the user. 

The user need not be concerned with array dimension statements. The program 
uses dynamic memory. Program dimensions and memory requirements are automatically 
allocated during execution. Moreover, they are automatically adjusted during a 



session to accommodate any other problem the user may choose to analyze. This 
capability is possible through the use of the Control Data Corporation (CDC) Common 
Memory Manager, which is a part of their operating system. 

The Common Memory Manager is a set of several FORTRAN-callable subroutines, each 
of which begins with the prefix CMM. These subroutines return to the calling program 
the first word address (FWA) of any array name. Another returns the offset needed 
to reposition this array at the end of the program, while a third routine sets aside 
that memory needed to execute the program. A fourth routine releases this space, and 
a fifth routine extends it. These routines are extremely useful within a variable 
dimension program like IMAGES and are used extensively therein. 

All array variables are listed in table III along with their memory assign- 
ments. This table provides the user with a quick index to each program variable and 
its position within the memory table. The variable name A is used extensively 
within both IMAGES and its executive subroutine MODAL. Other variables are listed in 
table III along with their algebraic equivalents. These variables can be correlated 
with those algebraic expressions discussed earlier. 

The "computational mode" is much more complex than either of the other modes. 

It is the backbone of IMAGES. All computations are initiated from this module. Some 
options require the user to choose between other suboptions, as can be seen in fig- 
ure 5. In some instances, the user is required to provide additional input. For 
example, the synthesis option will usually require the user to provide both the modal 
and the A matrices. In fact, the spectral synthesis returns to these inputs 
repeatedly until the specified tolerance a is satisfied. The "response" option is 
another example where additional input is required. In this instance, the user must 
specify either a step input or initial condition input as well as other parameters 
such as time increments. Other options will simply provide the specified computa- 
tional result. The user should be cautioned, however, to save needed results in all 
instances before entering the response option. In fact, the user should save results 
after each operation if it is desired to retain the data from completed steps for 
future use. 

The response option is provided as an analysis aid only. This option will pro- 
vide the user with a graphical display of system response to various types of input. 
These results may not be saved by the user. They are intended only as design aids. 
Hence, a particular design can be modified repeatedly until the desired response is 
obtained. 

The program is designed to provide the user with considerable numerical results. 
Hence, automatic scrolling could present difficulties for some terminals. Line limi- 
tations vary considerably from terminal to terminal. This is controlled internally, 
to some extent, within the program by using a PAUSE command. This command stops 
execution until the user enters a carriage return. The user can employ this command 
to control the rate of execution. It also allows the user to examine the returning 
output. If the numerical results are not progressing as expected, the user can 
return to an earlier step to modify his input. 

An example session is given in the appendix. It is not an example of a complete 
design problem. This example represents the end result of several sessions where the 
desired eigenvalues and eigenvectors have been refined to achieve the desired 
results. The user is encouraged to try this example. The program is extremely fast 
and requires less than 70K octal memory words to run. Response time is essentially 
instantaneous for most options. The various diagnostics will quickly become second 


13 



nature with only a few practice sessions such as this. Moreover, it will provide a 
valuable check on program compatibility. The reader is referred to the references 
for that logic used in selecting desired eigenvalues and eigenvectors in a practical 
design problem. 


CONCLUDING REMARKS 

An interactive digital computer program for modal analysis and gain estimation 
through eigensystem synthesis has been described. The program is modular in design 
and is extremely flexible. It is capable of both modal and spectral synthesis of 
multi-input control systems. The numerous options along with the incorporation of 
both memory management and random access have been combined to produce a truly user 
friendly, scratchpad program for the design of multi-input control systems. 

The program is very fast and requires less than 70K octal memory words to run. 
Some facilities impose limitations on programs to maximize the number of jobs which 
can be processed. Hence, these characteristics can be an added advantage to poten- 
tial users. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 30, 1983 
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APPENDIX 


EXAMPLE OF AN INTERACTIVE SESSION 


13.43.21 . RETURN (IMAGES) 

13.43 .22 . ATTACH (FTNMLIB/UN=LIBRARY) 

13 .43 .26 . REWIND (LMAP) 

13.43.27. LDSET(LIB=FTNMLIB,PRESETA=INDEF,MAP=SBEXXLMAP) 


MASTER MENU 
CHOOSE OPTION: 

1— BEGIN PROBLEM 

2— ENTER INPUT 

3— ENTER COMPUTATIONAL MODE 

4— SAVE PROBLEM 

5— RENAME FILE 

6— DUPLICATE FILE 

0— STOP 


? 1 

THERE ARE 2 PROBLEMS ON FILE. THESE ARE AS FOLLOWS: 
SAMPLE TEST 


IS THIS A NEW EXAMPLE? 
CHOOSE ONE: 

1— NEW 

2— OLD 

3— NEITHER 

? 1 


ENTER SEVEN (7) OR LESS CHARACTER FILE NAME: 
? sample 


PROBLEM NAME IS ON FILE!! 
DO YOU WISH TO PURGE? 
CHOOSE ONE:? 

1— YES 

2— NO 


? 2 


IS THIS A NEW EXAMPLE? 
CHOOSE ONE: 

1— NEW 

2— OLD 

3— NEITHER 

? 1 
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ENTER SEVEN (7) OR LESS CHARACTER FILE NAME: 
? example 


MASTER MENU 
CHOOSE OPTION: 

1— BEGIN PROBLEM 

2— ENTER INPUT 

3— ENTER COMPUTATIONAL MODE 

4— SAVE PROBLEM 

5— RENAME FILE 

6— DUPLICATE FILE 

0— STOP 


? 2 


PROBLEM IS A NEW ONE. ARRAYS MUST BE DIMENSIONED! 

DO YOU WANT A GRAPHIC DISPLAY OF MATRICES & THEIR DIMENSIONS 
CHOOSE ONE: 

1— YES 

2— NO 


? 1 

*** 

***** 

*** 


* 

* 


* * * 


***** 

* 

* 

* 

* 


* * 


* 

* 

* 

* 

* 

* 


* * * * 

* 

* 

* 

* 

* 

* 

* 

****** 

* * 

* 

* 

* 

* 

* 

* 


* * * 

***** 

**** 

* 

* 

* * 

****** 

***** * * 

* 

* 

* 

* 

* 

* 

* 


* * * * 

* 

* 

* 

* 

* 

* 

* 


* * 


* 

* 

** 

** 

* 

* 


* * * * 


***** 

***** 


N X 1 

N X N 

N X 1 

N X M 

M X 1 


SYSTEM 

STATE 

INPUT 

INPUT 

VECTOR 

MATRIX 

VECTOR 

MATRIX 

VECTOR 

PLEASE ENTER 

THE NUMBER OF 

STATE ELEMENTS, N. 
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THE CONTROL LAW MAY BE EITHER AN INPUT OR AN OUTPUT 
CONTROL: 


* * 

* 

* 

* 

* 


* * 


* * 

* 

* 

* 

* 


* * 


* * 

* 

* 

* 

* 

* 

* * 


* * 

****** * 

* 

* 

* 

* 

* * 


* * 

* 

* 

* 


***** 

* * 


* * 

****** ** 

* 

* 

* 

* 

* * 


* * 

* 

* 

* 

* 

* 

* * 

**** 

** ** 

* 

* 

* 

* 


** *« 

* 

***** 

* 

* 

* 

* 


***** 

* 








**** 








* 








* 

M X 1 


M X N 

N 

X 1 


M X 1 

* 

INPUT 

FEEDBACK 

STATE 


INPUT 


VECTOR 


MATRIX 

VECTOR 


VECTOR 


OR 









* 

* 

* * 

* * 


* 

* 


* 

* 

* * 

* * 


* 

* 


* 

* 

* * 

* * 

* 

* 

* 


* 

* 

****** * * 

* * 

* 

* 

* 


* 

* 

* * 

* 

***** 

* 

* 


* 

* 

****** ** * 

* 

* 

* 

* 


* 

* 

* * 

* 

* 

* 

* 

**** 

** ** 

* * 

* 


** ** 

* 

***** 

* * 

* 


***** 

* 








**** 








* 








* 


M X 1 

M X (R =N)(R =N) X 1 



M X 1 

* 


INPUT 

FEEDBACK 

OUTPUT 



INPUT 



VECTOR 

MATRIX 

VECTOR 



VECTOR 



BUT, THE NUMBER OF INPUTS WILL REMAIN UNCHANGED. 
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PLEASE ENTER THE NUMBER OF INPUTS, M. 


? 2 


PLEASE ENTER THE NUMBER OF OUTPUTS, R. 


? 4 


DESCRIBE YOUR INPUT (CHOOSE ONE) 

1— THE SYSTEM MATRICES A & B 

2— OUTPUT MATRICES (I.E., Y=CX+DU) 

3— BOTH "1" &"2" 

4 — EIGENVALUE/EIGENVECTOR MATRICES 

5— MODAL COUPLING MATRIX 

6— ENTER COMPUTATIONAL MODE 

7— SAVE INPUT 

8— PRINT A & B MATRICES 

9— PRINT C & D MATRICES 
0— EXIT 


? 1 


. . .EIGENVALUE/EIGENVECTOR ASSIGNMENT ALGORITHM 


...SYSTEM PARAMETERS... 
NO OF STATES = 4 

NO OF INPUTS = 2 

READ SYSTEM MATRICES A AND B 
ENTER A MATRIX 
ENTER BY ROWS: 


? -.56758, .859871 ,-2.192935,0.0 
? -.006922, -.147531 ,.230069,0.0 
? • 195034, -.970592,-. 079899, . 123380 
? 1 .0, .140541 ,0.0, 0.0 
ENTER B MATRIX 
ENTER BY ROWS: 
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7 . 829520 , .210565 
? -. 013460, -.115903 
? -. 006502, .0110501 
? 0 . 0 , 0.0 


: A : MATRIX 


-.567580 

.859871 

- 2.192935 

0.000000 

-.006922 

-.147531 

.230069 

0.000000 

.195034 

-.970592 

-.079899 

.123380 

1 .000000 

.140541 

0.000000 

0.000000 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

7 


: B : MATRIX 


.829520 

.210565 

-.013460 

-.115903 

-.006502 

.011050 

0.000000 

0.000000 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

9 
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DO YOU LIKE A & B? 
CHOOSE ONE: 

1— YES 

2— NO 


DO YOU WISH TO PERMUT A & B MATRICES? 
CHOOSE ONE: 

0— NO 

1— YES 


? 0 


DESCRIBE YOUR INPUT (CHOOSE ONE) 

1— THE SYSTEM MATRICES A & B 

2— OUTPUT MATRICES (I.E., Y=CX+DU) 

3— BOTH "1" &"2" 

4 — EIGENVALUE/EIGENVECTOR MATRICES 

5— MODAL COUPLING MATRIX 

6— ENTER COMPUTATIONAL MODE 

7— SAVE INPUT 

8 — PRINT A & B MATRICES 

9— PRINT C & D MATRICES 

0— EXIT 


? 7 


DESCRIBE YOUR INPUT (CHOOSE ONE) 

1— THE SYSTEM MATRICES A & B 

2— OUTPUT MATRICES (I.E., Y=CX+DU) 

3 — BOTH "1" &"2" 

4 — EIGEN VALUE/ EIGENVECTOR MATRICES 

5— MODAL COUPLING MATRIX 

6— ENTER COMPUTATIONAL MODE 

7 — SAVE INPUT 

8— PRINT A & B MATRICES 

9— PRINT C & D MATRICES 
0— EXIT 


? 0 


MASTER MENU 
CHOOSE OPTION: 

1— BEGIN PROBLEM 

2— ENTER INPUT 

3— ENTER COMPUTATIONAL MODE 

4— SAVE PROBLEM 

5— RENAME FILE 

6— DUPLICATE FILE 
0— STOP 
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COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 1 

CURRENT EPS IS 0. WHAT IS ITS NEW VALUE? 

? 0.000001 


: B1 INVERSE : 


1.242133 2.256626 

-.144251 -8.889970 

PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

7 

:DETERMINANT OF B1 = - .933 1 OE-O 1 

: S : MATRIX 

-.009670 -.112908 

0.000000 0.000000 

PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

7 
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: F : MATRIX 


.075129 -.123380 

0.000000 0.000000 

PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION 
PAUSE 

7 

: G : MATRIX 


.188764 -.978934 


1.000000 .140541 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION 
PAUSE 

7 


COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 5 
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COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 2 


CHOOSE ONE: 

1— SPECTRAL SYNTHESIS 

2— MODAL SYNTHESIS 

3— EXIT 


? 2 

CHOOSE ONE: 

0— MODAL MATRIX 

1— MODAL COUPLING MATRIX 


? 0 

SYSTEM MODES (AND THE OPTIONAL MODAL MATRIX) ARE REQUIRED. 
ENTER DESIRED MODES 


A MAXIMUM OF 4 REAL MODES ARE EXPECTED 
COMPLEX MODES WILL AUTOMATICALLY DIMINISH THIS NUMBER 


ENTER EACH IN TABLE FORM WITH COMMA (,) SEPARATORS: 


TYPE , REAL PART , IMAG PART 

SPACING IS NOT IMPORTANT AND TYPE CODES ARE AS FOLLOWS: 

0 — IMPLIES REAL MODE; 1 — IMPLIES COMPLEX MODE; 2 — ENDS INPUT 


? 0 ,-. 72 , 0.0 
? 1 ,-.15, .4 
? 0 ,-. 031 , 0.0 
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DESIRED MODES 

TYPE REAL PART IMAG PART 

0 -.7200E+00 0. 

1 1500E+00 .4000E+00 

0 - .3 100E-01 0. 

CHOOSE ONE: 

1— CORRECT MODES 

2— CONTINUE 


? 2 


ENTER THE M X N DESIGN MATRIX " Z " BY ROWS 


DO YOU WANT A GRAPHIC DISPLAY OF MATRICES & THEIR DIMENSIONS 
CHOOSE ONE: 

1— YES 

2— NO 


? 1 


***** 

* 

* * * 

* 

* * * 
****** * 

* * * 
****** * 

* * * 

* 

* * 

* 

***** 


********* 

* 


* 


* 


* 

* 

* 


* 

********* 


** 

** 


** 

** 


***** 

* 

* * * 
* 
* 

* * * 

* * 

* 

* * * * * 

* 

* * * 
* 

***** 


N X 1 
EIGEN- 
VECTOR 
ENTER ROW 1 
? 1 .0,0. 0,0.0, -.045841 
ENTER ROW 2 
? 0.0, .5, .5, .10927 


M X 1 
SPECIFIED 
ELEMENTS 


(N-M) X 1 

ASSIGNED 

ELEMENTS 


***** 

* 

* 

* 

* 

* 
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DESIGN PARAMETER MATRIX " Z " 


1 .000000 
0.000000 


0.000000 

.500000 


0.000000 

.500000 


-.045841 

.109270 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

? 

CHOOSE ONE: 

1— CORRECT "Z" 

2— CONTINUE 


COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 5 


COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 3 
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FEEDBACK GAINS 
CHOOSE ONE: 

1— ENTER STATE FEED BACK GAINS 

2— CALCULATE FEEDBACK GAINS 

3— EXIT 


CHOOSE SYSTEM: 

1— STATE FEEDBACK 

2— OUTPUT FEEDBACK 


MATRIX OF EIGENVALUES 


-.720000 

0.000000 

0.000000 

0.000000 

0.000000 

-.150000 

.400000 

0.000000 

0.000000 

-.400000 

-.150000 

0.000000 

0.000000 

0.000000 

0.000000 

-.031000 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION 
PAUSE 

7 


: V : MATRIX (MODAL MATRIX) 


.584162 

0.000000 

0.000000 

-.045841 

0.000000 

.269093 

.269093 

.109271 

-.022072 

-.579930 

.709355 

.137647 

-.811337 

.051806 

-.113974 

.983368 
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PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION. 

PAUSE 

7 


DETERMINANT OF : V : = 

.19181E+00 



: V-INVERSE : 




1 .827418 

.007014 

.010662 

.082917 

-.203078 

2.034071 

-.791655 

-.124680 

-.394926 

1 .645819 

.736485 

-.304382 

1 .472654 

.089380 

.135863 

1 .056615 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

? 

NUMERICAL CONDITION OF SOLUTION 
(NORM OF RECIPROCAL VECTORS) 

MODE CONDITION NUMBER 


-.72000E+00C+/-) 

J 

0 . 

.18293E+01 

-. 15000E+00(+/-) 

J 

.40000E+00 

.28846E+01 

-.31000E-01(+/-) 

J 

0 . 

.18198E+01 

+ BK : MATRIX 





-.766513 

-.002823 

-.004291 

-.033373 

-.001501 

-.190628 

.166251 

-.005603 

.196346 

-.957383 

-.093858 

.124335 

1 .000000 

.140541 

.000000 

.000000 
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PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED .PRESS 'ENTER' TO CONTINUE EXECUTION. 
PAUSE 

9 


: K : MATRIX 


-.234868 -1.168835 2.574575 -.054098 

-.019499 .507576 .251622 .054626 


PROGRAM EXECUTION IS TEMPORARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION 
PAUSE 

7 

: ERROR MATRIX : 


-.000000 .000000 .000000 -.000000 

-.000000 0.000000 .000000 .000000 


PR °PAUSE EXECUTI0N IS TEMP0RARILY INTERRUPTED. PRESS 'ENTER' TO CONTINUE EXECUTION. 


...MESSAGE ... 


IF ELEMENTS OF ERROR MATRIX HAVE NOT APPROACHED ZERO 
SIGNIFICANT ROUND OFF ERRORS HAVE OCCURRED IN COMPUTATION. 

PROBABLE CAUSES ...1) MATRIX B1 IS ILL CONDITIONED 

TRY PERMUTING SYSTEM MATRICES TO GENERATE BETTER B1 MATRIX 

2) MODAL MATRIX IS ILL CONDITIONED 
CHANGE DESIGN VECTORS 
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COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— ' TIME RESPONSE 

5— SAVE 

6— EXIT 


? 5 


COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3— GAINS 

4— ' TIME RESPONSE 

5— SAVE 

6— EXIT 


? 4 


CHOOSE OPTION: 

0— EXIT 

1— STATE RESPONSE 

2— STATE+CONTROL RESPONSE 

3— STATE+OUTPUT RESPONSE 

4— STATE+OUTPUT+CONTROL RESPONSE 


? 1 


ENTER THE NUMBER OF TIME STEPS: 

? 50 

ENTER NUMBER OF OUTPUT RESPONSES, R 
? 4 

WHAT RUN NUMBER IS THIS? 

? 1 

ENTER THE TIME INCREMENT 

? .2 

CHOOSE THE DESIRED TYPE OF RESPONSE: 

0— INITIAL CONDITION RESPONSE 

1— STEP RESPONSE 
? 0 

ENTER THE " 4" X " 1" INITIAL STATE VECTOR 

? 0 . 0 . 0 . 0 ,-. 1 , 0.0 
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STATE VARIABLE RESPONSE 

TIME RESPONSE PARAMETERS 


NO OF STEPS = 50 

TYPE OF RESPONSE = 0 

TIME INCREMENT = .200 SECONDS 


...INITIAL CONDITION RESPONSE ... 


. . . I.C. VECTOR . . . 


0.000 0.000 -.100 0.000 
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CHART 1 


0.0000 3 


.4000 3 

.6000 3 

.8000 
1.0000 
1 .2000 

1 .4000 
1.6000 
1.8000 
2.0000 
2.2000 

2.4000 
2.6000 
2.8000 
3.0000 
3.2000 

3.4000 
3.6000 
3.8000 

4.0000 

4.2000 

4.4000 

4.6000 

4.8000 

5.0000 

5.2000 

5.4000 

5.6000 

5.8000 
6.0000 
6.2000 

6.4000 
6.6000 
6.8000 
7.0000 
7.2000 

7.4000 
7.6000 
7.8000 
8.0000 
8.2000 

8.4000 
8.6000 
8.8000 
9.0000 
9.2000 

9.4000 
9.6000 
9.8000 


3 


2 4 
41 
41 
41 
41 
41 
41 
1 


2 

2 


2 

2 

2 

2 
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2 

2 

2 


4 

4 

4 


2 

2 

2 


4 

4 

4 

4 

4 

43 


2 

2 


4 

4 

4 

4 

4 

4 

4 

> 4 
4 

42 

42 

4 2 

4 ; 

4 ; 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

3 

T 3 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

! 1 
21 
21 
2 
2 
12 
12 
12 
1 2 
1 2 
1 2 


3 

3 

3 


3 

3 

3 

3 



-.1000 


-.0815 


-.0630 -.0445 


-.0260 


-.0075 


.01 10 


.0295 
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CHOOSE OPTION: 

0— EXIT 

1 — STATE RESPONSE 

2— STATE+CONTROL RESPONSE 

3— STATE+OUTPUT RESPONSE 

4 — STATE+OUTPUT+CONTROL RESPONSE 


? 0 


COMPUTATIONAL MODE: 

CHOOSE ONE— 

1— CONSTRAINING RELATION 

2— SYNTHESIS 

3 — GAINS 

4— TIME RESPONSE 

5— SAVE 

6— EXIT 


? 6 


MASTER MENU 
CHOOSE OPTION: 

1— BEGIN PROBLEM 

2— ENTER INPUT 

3— ENTER COMPUTATIONAL MODE 

4— SAVE PROBLEM 

5— RENAME FILE 

6— DUPLICATE FILE 
0— STOP 


? 0 

HAS THE FILE BEEN SAVED? 
(YES/NO) 

? yes 

14.16.55 . REWIND(LGO) 

14.16.56. RETURN (FTNMLIB) 

14. 16. 57. REPLACE, STRTUP. 
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TABLE I.- CONTROL WORDS AND POSSIBLE VALUES 


Name 

Function or description 

Possible values 

EPS 

Tolerance parameter 

If |det(A) | < EPS, the matrix A is assumed to be 
singular 

IENTRY 

Control parameter 

-2 - Computes time response for a specified V 
and A, 

-1 - Computes eigenvalues/eigenvectors of matrix A, 
numerical condition (sensitivity) of modal 
matrix, and time response 

0 - Normal setting for both modal and spectral 

synthesis 

1 - Normal setting for gain estimation 

2 - Setting for feedback synthesis 

INPUT 

Controls type of response 

0 - Initial condition response 

1 - Step response 

IPRMUT 

Reorders state variables to 

0 - No 


improve numerical 
condition of B 1 matrix 

1 - Yes 

IPUNCH 

Controls print 

0 - No 

1 - Yes 

ITSOLN 

Controls time response 

0 - No time response 

1 - State variable response 

2 - State + control response 

3 - State + output response 

4 - State + control + output response 

MODE(I) 

Type of mode (complex modes 

0 - Real mode 


occur in pairs) 

1 - Complex mode 

NPHASE 

Control parameter 

-2 - Permutes A & B and C & D 

-1 - Computes B^ 1 , F, S, and G matrices 

0 - Selects modal synthesis 

1 - Computes modal coupling matrices C^ for each 

specified mode 

2 - Selects spectral synthesis 

NRUN 

Identification number 
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TABLE II.- DESCRIPTION OF SUBROUTINES 


Name 

Function or description 

CONDN 

Computes numerical ill-conditioning (sensitivity) of modal 
matrix V 

DIAG 

System routine for generating the quasi-diagonal equivalent 
of the eigenvalue and eigenvector matrices 

EIG 

1 

System routine for finding eigenvalues and eigenvectors i 

ENT 

Reads a matrix j 

LUDECP & SUBST 

Solves a system of linear equations 

MATEQN 

Inverts a square matrix 

MINV 

j 

System routine for inverting matrices j 

j MULT 

i 

Multiplies two matrices 

NORM 

Normalizes a vector 

OGAIN 

i 

Computes the output feedback gains approximation (ic) to a state' 
feedback design (K) j 

PERMUT 

j 

Permutes rows and columns of a matrix | 

l 

PLOT 

! 

System routine for generating printed plots 

| PRT 

Prints a matrix 

SGAIN 

Computes the equivalent state feedback gains (K) for an output ! 
feedback design i 

j SPECTRL 

j 

Interactive subroutine for eigensystem placement 

! SYNTHS 

Computes modal coupling matrices and eigenvectors for specified 
\ and z, (i = 1, 2, n) 

TFUNC 

Evaluates e^ 

TRANSF 

Computes the closed-loop system after applying output feedback 

USOLN 

Computes actuator response u = K + u 

x ref 

XSOLN 

Computes state variable response time histories for 
x ( t ) = Ve At V~ 1 x(0) 

YSOLN 

Computes designated output response y = C x + D u 

o o o 
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Variable 

Algebraic 

equivalent 

Memory assignment 

First word 
address 

Length 

Upper 

limit 

A 

A 

IOFF 

N 

X 

N 

11 

AE 


11 3 

N 

X 

N 

114 

AHAT 

A 

11 2 

N 

X 

N 

11 3 

B 

B 

11 

N 

X 

M 

12 

B1 

B, 

12 

M 

X 

M 

13 

B2 

B 2 

13 

L 

X 

M 

14 

BINV 

s ? 1 

14 

M 

X 

M 

15 

BTEMP 


15 

N 

X 

M 

16 

C 

C 

117 

N 

X 

N 

11 8 

CCT 

Im{C. } 

147 

LL 

X 

MM 

148 

CH 


157 

K 

X 

N 

158 

Cl 


159 

K 

X 

K 

160 

CM 

Im{\^ } 

136 

NM 


137 

CMK 

Im{M (k) } 

178 

KK 

X 

KK 

179 

CO 


150 

K 

X 

N 

151 

CP 


14 2 

LL 

X 

LL 

143 

CQ 


143 

LL 

X 

MM 

144 

CT 


158 

K 

X 

K 

159 

CW 


16 2 

N 

X 

K 

163 

CY 


169 

NO 

X 

N 

170 

D 

D 

19 

N 

X 

N 

11 0 

DG 


155 

K 

X 

N 

156 

DO 


151 

K 

X 

M 

152 

DT 


156 

K 

X 

M 

157 

E 


114 

L 

X 

N 

11 5 

El 


17 5 


N 


176 

ER 


174 


N 


175 

F 


115 

N 

X 

N 

116 

FC 


185 

2 

X 

M 

186 

FR 


19 2 


N 


19 3 

FRC 

r 

188 

2 

X 

NM 

189 

G 

K 

11 8 

M 

X 

N 

119 

G1 

K 1 

119 

M 

X 

M 

120 

G2 

1 

120 

M 

X 

L 

121 

GH 


15 2 

M 

X 

M 

15 3 

GO 


15 3 

M 

X 

K 

154 

GOT 


154 

M 

X 

K 

155 

GRT 

[g (k) ]T 

181 


L 


184 

GRTC 


18 7 

2 

X 

KK 

188 

GT 


121 

M 

X 

N 

122 










TABLE III 


Continued 


Variable 

Algebraic 

equivalent 

Memory assignment 

First word 
address 

Length 

Upper 

limit 

H 


116 

N 

X 

N 

117 

HC 

h (k-1) 

r 

19 3 

2 

X 

NM 

194 

HR 

18 3 


NM 


184 

HRC 


189 

2 

X 

NM 

190 

IK 


160 


K 


161 

IL 


141 


L 


142 

ILL 


148 


LL 


149 

IM 


124 


M 


125 

IN 


122 


N 


123 

INO 


171 


NO 


172 

IRCCK 


184 


2L 


185 

IRCK 


182 


L 


183 

JK 


161 


K 


16 2 

JLL 


149 


LL 


150 

JM 


125 


N 


126 

JN 


123 


N 


124 

MODE 


134 


NM 


135 

Q 


179 

L 

X 

L 

180 

QC 

g(k-1) 

180 

KK 

X 

KK 

181 

QR 

191 

N 

X 

N 

19 2 

R 

F 

16 

L 

X 

L 

17 

RE 

Re{X i } 

135 


NM 


136 

RM 

Re{M^ k ) } 

19 0 

N 

X 

N 

191 

RMK 

177 

L 

X 

L 

178 

RRT 

Re{Q„} 

140 

L 

X 

M 

141 

S 

S 

17 

L 

X 

M 

18 

SA 

A 11 

130 

M 

X 

M 

131 

SB 

A 21 

131 

M 

X 

L 

13 2 

SC 

132 

L 

X 

M 

133 

SD 


133 

L 

X 

L 

134 

T 

G 

18 

L 

X 

M 

19 

U 


11 0 

N 

X 

N 

11 1 

UC 


165 


M 


166 

UO 


164 


M 


165 

UPLT 


167 


LU 


168 

UREF 


172 


M 


17 3 

V 


11 1 

N 

X 

N 

11 2 

VC 


186 


NN 


18 7 

WC 

Im^} 

146 

MM 


147 

WK 

Re {w i } 

176 


NN 


177 

WR 

138 


L 


139 









Memory assignment 


Variable 


X 

XI 
XO 

XPLT 

XR 

YC 

YO 

YPLT 

YR 

YU 

Z 

ZC 

ZR 


Algebraic 

equivalent 

First word 
address 

Length 

Upper 

limit 


126 


127 


129 


130 


127 


128 


166 


167 


128 


129 


144 

11 ^ • •' 

145 


173 


174 


168 


169 


137 


138 


170 

NO 

171 


163 

N 

164 

Im{z. } 

146 

MM 

147 

Re{z i } 

139 

M 

140 












Figure 1.- Flow chart for input module. 

















Figure 2.- Computational flow for eigensystem synthesis. 

















Figure 4.- Flow of response computations. 
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OLD 
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1-YES 

1-YES 



2-NO 

2— NO 


t OLD/NEW 
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STATE ELEMENTS. N. 
INPUTS. M. 

OUTPUTS. R. 


[DESCRIBE YOUR INPUT (CHOOSE ONE) 

1 - THE SYSTEM MATRICES A & B 

2- OUTPUT MATRICES (I.E. Y=CX+DU) 

3- BOTH "1" & "2" 

.4— EIGENVALUE/EIGENVECTOR MATRICES 

5- MODAL COUPLING MATRIX 

6- ENTER COMPUTATIONAL MODE 

7- SAVE INPUT 

8- PRINT A & B MATRICES 
PRINT C & D MATRICES 

0-EXIT 


NEW 


NEITHER 


PURGE/NEW 


I ENTER SEVEN (7) OR LESS CHARACTER FILE NAME-1 

_L 
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| IS THIS A NEW EXAMPLE? 
CHOOSE ONE: 

1— NEW 

2— OLD 

3— NEITHER 


OLD 


PROBLEM NAME IS NOT ON FILE 1 ' 
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EXI T I I.E. , 0)[ MODIFY 

| PROBLEM IS AN OLD ONE. CHOOSE ONE- 

1— MODIFY INPUT 

2- RETURN 
i RETURN 


-PROBLEMS ON FILE. THESE ARE AS FOLLOWS -| 


|MASTER MENU 
CHOOSE OPTION: 

1- BEGIN PROBLEM 

2- ENTER INPUT 

3- ENTER COMPUTATIONAL MODE] 

4- SAVE PROBLEM 

5- RENAME FILE 

6- DUPLICATE FILE 
O-STOP 


© 


© 


HAS THE FILE BEEN SAVED? 
.(YES /NO) 


1 1 S THIS AN OVERSIGHT? < YES/NO) 
J 


CURRENT FILE NAME IS 
WHAT IS THE NEW NAME? 

ENTER SEVEN (7 ) OR LESS CHARACTER FILE NAME: 


STOP 


1 OR 3 


[READ SYSTEM MATRICES A AND B 
ENTER A MATRIX 
[ENTE R BY ROWS : 

ENTER B MATRIXl 
ENTER B ROWS: 


CORRECTIONS FOR: A: MATRIX [ 
ENTER ROW & COLUMN (RR. CC ): I 


ENTER NEW ELEMENT | - 


|DOYOU LIKE A&B ?| 
CHOOSE ONE: 

1- YES 

2- NO 



T 1 OR 3 



3 


WHICH MATRIX? 
CHOOSE ONE: 

1- A MATRIX 

2— B MATRIX 





2 

CORRECTIONS FOR : B: MATRIX 

2 


ENTER ROW & COLUMN! RR CCF 


3-BOTH 





IT 


DO YOU WISH TO PERMUTE A & B MATRICES? 
CHOOSE ONE: 


| COLUMN PERMUTATION} ^- 


0- NO 

1- YES 


[WHICH MATRIX? 
CHOOSE ONE: 

1— A MATRIX 

2— B MATRIX 

3— BOTH 


ENTER PERMUTATION SEQUENCE FOR : B: MATRIX 
ROW PERMUTATION 

' n 


10R 3 


I COLUMN PERMUTATIONl - 


1 


ENTER PERMUTATION SEQUENCE FOR- A: MATRIX 
ROW PERMUTATION 


©j PLEASE ENTER THE NUMBER OF OUTPUTS - ^ 


ENTER C MATRIX 
ENTER BY ROWS : 


ENTER D MATRIX 
ENTER BY ROWS: 


I 


DO YOU LIKE C & D ?- 
CHOOSE ONE: 

1- YES 

2- NO 
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CHOOSE ONE 
O-NO 
1-YES 


'CORRECTIONS FOR: C: MATRIX 

ENTER ROW & COLUMN (R R. CC ):| 
I 1 OR 3 
WHICH MATRIX? 

CHOOSE ONE: 

1— C MATRIXl 

2— D MATRIXl 

3— BOTH 


] ENTER NEW ELEMENT \ 


CORRECTIONS FOR : D: MATRIX 
ENTER ROWS COLUMN (RR. CC): 


T 


I ENTER NEW ELEMENT I 


I COLUMN PERMUTATION^ 


[WHICH MATRIX? 
CHOOSE ONE: 

1— C MATRIX 

2— D MATRIX 

3— BOTH 


1 OR 3 


-ENTER PERMUTATION SEQUENCE FOR: D- MATRIX 
ROW PERMUTATION 

r , M 

| COLUMN PERMUTATIONl - 

ENTER PERMUTATION SEQUENCEFOR: C : MATRIX) 
ROW PERMUTATION 


Figure 5.- Flow chart of iterative discourse. 
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Figure 5.- Concluded. 
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